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Abstract 

We present the results of a realistic global fit of the Lagrangian parameters of the Min- 
imal Supersymmetric Standard Model assuming universality for the first and second 
generation and real parameters. No assumptions on the SUSY breaking mechanism are 
made. The fit is performed using the precision of future mass measurements of super- 
partners at the LHC and mass and polarized topological cross-section measurements 
at the ILC. Higher order radiative corrections are accounted for whereever possible 
to date. Results are obtained for a modified SPSla MSSM benchmark scenario but 
they were checked not to depend critically on this assumption. Exploiting a simulated 
annealing algorithm, a stable result is obtained without any a priori assumptions on 
the values of the fit parameters. Most of the Lagrangian parameters can be extracted 



at the percent level or better if theoretical uncertainties are neglected. Neither LHC 
nor ILC measurements alone will be sufficient to obtain a stable result. The effects of 
theoretical uncertainties arising from unknown higher-order corrections and parametric 
uncertainties are examined qualitatively. They appear to be relevant and the result 
motivates further precision calculations. The obtained parameters at the electroweak 
scale are used for a fit of the parameters at high energy scales within the bottom-up 
approach. In this way regularities at these scales are explored and the underlying 
model can be determined with hardly any theoretical bias. Fits of high-scale parame- 
ters to combined LHC+ILC measurements within the mSUGRA framework reveal that 
even tiny distortions in the low-energy mass spectrum already lead to inacceptable \ 2 
values. This does not hold for "LHC only" inputs. 



1 Introduction 

Provided low-energy Supersymmetry (SUSY) [1] is realized in Nature, the next generation of 
colliders, the Large Hadron Collider (LHC) [2] and the International Linear Collider (ILC) [3] 
are likely to copiously produce SUSY particles and will allow for precise measurements of 
their properties. Once SUSY is established experimentally, it is the main task to explore 
the unknown mechanism of SUSY breaking (SSB). While specific models of SUSY breaking, 
like e. g. minimal supergravity (mSUGRA) [4] can be tested against the data in a relative 
straight-forward manner (see e. g. [5]), an exploration of the parameters of the general 
Minimal Supersymmetric Standard Model (MSSM) [6] parameter space is significantly more 
ambitious. A first step in this direction has been done in [7] where a fit of the real and 
flavour-diagonal MSSM parameters has been performed using MINUIT [8]. The use of 
MINUIT implied some restrictions, e. g. the starting point for the fit has to be close to 
the real values and also the correlation matrix depends significantly on the starting point 
due to the large number of parameters. These limitations can be overcome by the methods 
presented in [9] and it is among the aims of this paper to demonstrate this in detail. 

The general MSSM Lagrangian contains more than 100 new parameters which are only 
related to each other through the unknown SSB mechanism. Many of these new parame- 
ters are CP-violating complex phases which are to some extent limited by the absence of 
neutron and electron electric dipole moments [10] and flavour-nondiagonal couplings which 
are bounded by the absence of flavour-changing neutral currents both in the quark and lep- 
ton sector [11, 12]. It seems therefore justified to initially consider a somewhat constrained 
MSSM with real parameters and flavour- diagonal couplings. Furthermore universality of the 
first and second generation parameters appears to be a reasonable approximation while large 
Yukawa couplings lead to significant differences for the third generation. Applying these con- 
straints, the number of free parameters is reduced to 19, including the top quark mass as a 
parameter to account for parametric uncertainties. While still many, these parameters may 
be confronted with an even larger number of independent observables at the LHC and ILC. 
Thus it is of high interest how well this 19-dimensional parameter space can be restricted 
with future measurements. 
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The experimental collabarations ATLAS [13] and CMS [14] at the LHC have performed 
detailed simulations of the possibilities to extract mass information from their future data, 
predominantly from the reconstruction of kinematic end-points of involving leptons and jets 
in cascade decays of gluinos and squarks [15] produced in proton proton collisions at 14 TeV 
centre-of-mass energy. 

At a future electron positron collider for collision energies up to 1 TeV, the International 
Linear Collider (ILC), the kinematically accessible part of the superpartner spectrum can be 
studied in great detail due to favourable background conditions and the well-known initial 
state [15]. 

The attempt of an evaluation of the MSSM Lagrangian parameters and their associated 
errors is only useful if the experimental errors of the future measurements are known and un- 
der control. The by far best-studied SUSY scenario to serve as a basis for such an evaluation 
is an mSUGRA-inspired benchmark scenario, the SPSla scenario. For this scenario with a 
relatively light superpartner spectrum, a wealth of experimental simulations exists and has 
recently been compiled in the framework of the international LHC/ILC study group [15]. 
In this paper we take a point close in parameter space which is consistent with low energy 
data and respects the dark matter constraints. We perform a global fit of the 19 param- 
eters to those expected measurements augmented by possible measurements of topological 
cross-sections (i. e. cross sections times branching fractions) that will be possible at the ILC 
with polarised beams. Since a detailed experimental simulation for some measurements is 
lacking, we estimate their uncertainties conservatively from the predicted cross-sections and 
transferring experimental efficiency from well-studied cases. 

While at leading order, certain subsets of the parameters of the Lagrangian only influence 
certain subsets of observables (e. g. only three parameters determine the masses of the 
charginos) this 'block-wise diagonal' mapping of parameters to observables no longer holds at 
loop-level. At loop-level, in principle each observable depends on each Lagrangian parameter, 
thus making the extraction of the parameters much more involved. This is particularly 
striking in the supersymmetric Higgs sector where due to the large third generation Yukawa 
couplings, radiative corrections to the mass of the lightest Higgs boson are in general of the 
order of 30 % to 50 %. However, also for the superpartner properties radiative corrections 
become important as soon as experimental precision enters the percent level which is clearly 
the case for many of the measurements possible at the LHC and in particular at the ILC. 

The calculation of higher-order corrections to SUSY observables has started in many sec- 
tors, see e. g. ref. [16] for an overview. However, a coherent framework for these corrections 
and a study of the transformation of the parameters defined in different renormalisation 
schemes has only started recently in the framework of the Supersymmetry Parameter Anal- 
ysis (SPA) project [17]. 

As theoretical basis for the global fit presented in this paper we use the calculations as 
implemented into the program SPheno [18]. In SPheno, the masses and decay branching 
fractions of the superpartners are calculated as well as production cross sections in e + e~ 
collisions. The calculation of the masses is carried out in the DR scheme and the formulae 
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for the 1-loop masses are used as given in ref. [19]. In the case of the Higgs boson masses, the 
2-loop corrections as given in ref. [20] are added. The calculation of the branching ratios is 
performed at tree-level using however running couplings evolved at the scale corresponding 
to the mass of the decaying particle. In case of the cross sections tree-level formulae are used 
except for the production of squarks in e + e~ annihilation where the formulae of ref. [21] 
are used. In addition we have added ISR corrections as given in [22]. For the evolution of 
parameters between various energy scales the 2-loop RGEs as given in [23] are used. 

In this paper we present the result of a global fit of the MSSM Lagrangian parameters at 
the electro-weak scale for a slightly modified SPSla benchmark scenario using the program 
Fittino [24, 9]. For a similar program see [25]. Previous evaluations of the errors of those 
parameters [7] did not attempt to develop a strategy to extract those parameters from data 
without a priori knowledge. Within Fittino, special attention is given precisely to this 
task, i. e. to find the parameter set which is most consistent with the data before a careful 
evaluation of errors and correlations is performed. The obtained parameters are then evolved 
to high energy scales, e. g. a GUT scale, taking into account all correlations between the 
errors. This is done within the bottom-up approach, e. g. no assumptions regarding the 
underlying high energy model are used. For further details see [7]. 

This paper is organized as follows. In Section 2 we shortly describe the approach used in 
Fittino. In Section 3, the modified SPSla benchmark scenario (SPSla') and assumptions for 
the input observables are explained. The results of the fit and the error evaluation method 
are summarized in Section 4. The extrapolation of the obtained fit parameters and their 
errors to high energy scales is carried out in Section 5. A fit within the mSUGRA framework 
is perfromed in Section 6 and conclusions are drawn in Section 7. 

2 Fit Procedure 

From the numerous fitting options provided by Fittino we have chosen the following fit proce- 
dure to extract the low-energy Lagrangian parameters. First start values for the parameters 
are calculated using tree- level relations between parameters and a few observables [9]. For 
fits with many parameters these values are not good enough to allow a fitting tool like 
MINUIT [8] to find the global minimum due to the amount of loop-level induced cross- 
dependencies between the individual sectors of the MSSM. Therefore, in a second step, the 
parameters are refined, using a simulated annealing approach [9, 26, 27]. As a result, the 
parameter values are close to the global minimum so that a global fit using MINUIT can 
find the exact minimum in a third step. To determine the parameter uncertainties and 
correlations many individual fits with input values randomly smeared around their true val- 
ues within their uncertainty range are carried out. The parameter uncertainties and the 
correlation matrix are derived from the spread of the fitted parameter values. 

The advantages of this approach are: 

• All available measurements can be used at once to extract the maximal possible infor- 
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mation from the data. 

• No a priori knowledge is required for the fit. The calculation of tree-level start values 
is a time-saving way to obtain reasonable start values for the fit. A brute force scan 
of the parameter space would not be feasible for a large set of parameters. 

• Correlations between input observables can easily be taken into account. 

The results presented in this article were obtained with Fittino version 1.1.1. More 
detailed information on the algorithms used by Fittino can be found in [9]. 



3 SPSla' Inspired Fit 

The SPSla' inspired scenario studied in this article is defined by the following high-scale 
parameters [17]: 

mo = 70 GeV universal scalar mass (1) 

m i/2 — 250 GeV universal gaugino mass (2) 

A = —300 GeV universal trilinear coupling (3) 

tan/5 = 10 (4) 

sign(/x) = +1 (5) 

As opposed to SPSla [28], SPSla' has the advantage that it is fully consistent with all 
available measurements including cosmological data [29]. 

Although in the definition of this scenario, gravity mediated SUSY breaking (mSUGRA) 
is assumed, no assumption on the SUSY breaking mechanism is made in the reconstruction 
of the Lagrangian parameters. This generality entails the introduction of many soft SUSY 
breaking parameters. The advantage of this approach is that the parameters are recon- 
structed at the low-energy scale without unnecessary assumptions and can subsequently be 
extrapolated to the high-scale to learn about the SUSY breaking mechanism ( "bottom-up" 
approach) [7]. 



3.1 Fit Assumptions 

The current version of Fittino (version 1.1.1) is able to fit all low-energy SUSY parameters 
of theories fulfilling the following properties: 

• There is no CP violation in the SUSY sector of the theory, i. e. all phases vanish. 

• No inter-generation mixing is present. 
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• Mixing within the first and second generation is zero. 

With these assumptions, 24 SUSY parameters remain from the initial 105 parameters in the 
general case. 

For the fit presented here, the 24 MSSM parameters available in Fittino have been further 
reduced. We assume that the first and second generation sparticle masses are almost degen- 
erate which is motivated by the fact that no deviation from SM predictions have been found 
up to now in low energy data, e. g. in kaon physics [11, 12]. In the squark sector, one unified 
squark mass parameter Mg L is assumed for the superpartners of the left-handed u, d, s and c 
quarks and one mass parameter Mg R is used for the superpartners of the right-handed light 
quarks. Using all assumtions, 18 free MSSM parameters are left. For the fit, the low-energy 
parameters calculated from equations 1 to 5 are slightly modified so that this unification is 
exact. Additionally, the top quark mass m t is fitted to account for parametric uncertainties. 

Instead of the trilinear couplings A t , Ay, and A T , the parameters 

X t = A t — fi/ tan/3 
X h = A h — fitanj3 
X T = A T — /it an/? 

are fitted in order to reduce correlations with tan (3 and \i. 

3.2 Input Observables to the Fit 

A number of anticipated LHC and ILC measurements serve as input observables to the fit. 
For the ILC, running at centre-of-mass energies of 400 GeV, 500 GeV and 1 TeV is considered 
with 80 % electron and 60 % positron polarization. The predicted values of the observables 
are calculated using the following prescriptions: 

• Masses: 

The experimental uncertainties of the mass measurements are taken from [15]. For 
the LHC, the mass information is extracted from measurements of edge positions in 
mass spectra. In the gaugino and Higgs sector the precision is driven by the ILC, for 
strongly interacting SUSY particles by the LHC. The benefit of combined analyses at 
LHC and ILC is taken into account. 

• Cross-sections: 

Only e + e~ cross-sections are included in the fit. However, the measurement of absolute 
cross-sections is impossible for many channels, in which only a fraction of the final 
states can be reconstructed. Therefore, absolute cross-section measurements are only 
used for the Higgs-strahlung production of the light Higgs boson, which is studied in 
detail in [30]. 
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(6) 
(7) 
(8) 



• Cross-sections times branching fractions: 

Since no comprehensive study of the precision of cross-section times branching fraction 
measurements is available, the uncertainty is assumed to be the error of a counting 
experiment with the following assumptions: 

— The selection efficiency amounts to 50 %. 

— 80 % polarization of the electron beam and 60 % polarization of the positron can 
be achieved. 

— 500 fb -1 per centre-of-mass energy and polarization is collected. 

— The relative precision is not allowed to be better than 1 % and the absolute 
accuracy is at most 0.1 fb to account for systematic uncertainties. 

All production processes and decays of SUSY particles and Higgs bosons are used 
which have a cross-section times branching ratio value of more than 1 fb in one of the 
e + e~ polarization states LL, RR, LR, RL at 500 GeV and LR or RL at y/s = 400 GeV 
and yfe = 1000 GeV. 

• Branching fractions: 

The four largest branching fractions of the lightest Higgs boson are included. The 
uncertainties on the Higgs branching fractions are taken from [30]. 

• Standard Model parameters: 

The present uncertainties for and mz are used as conservative estimates. The 
uncertainty for the top mass m t is assumed to be 50 MeV (experimental) [30] and 
100 MeV (theoretical) [31]. 

• Mixing angles: 

For tree-level estimates of tan/3, \i and the parameters of the gaugino sector the 
chargino mixing angles cos20j, and cos2</>/j are used. Their values are reconstructed 
using tree-level relations from chargino production cross-sections at different beam po- 
larizations [32, 33]. No use is made of those observables in the fit. The fit result is 
independent of their assumed uncertainties. 

In order to check the influence of theoretical uncertainties on the fit results, the fit has 
been performed twice, once with experimental uncertainties only and a second time with 
experimental and estimated theoretical uncertainties. To be conservative we take present 
theoretical uncertainties. As an estimate, the scale uncertainties as given in [17] are used for 
the mass predictions which also induce a shift in the 'edge' observables. They are obtained 
by varying the scale, where the parameters are decoupled from the RGE running and the 
shifts to the pole masses are calculated, from mz to 1 TeV. In this way one gets an estimate 
of the missing higher order corrections for the shift from the running mass to the pole 
mass. Here we use the complete 1-loop formulae given in [19] for the SUSY masses and in 
addition the 2-loop corrections for the Higgs sector as given in [20]. The experimental and 
theoretical contributions are added in quadrature. The assumed experimental and theoretical 
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Figure 1: Relative uncertainties of the parameter measurements with (red, left error bars) 
and without (blue, right error bars) theoretical uncertainties. 

uncertainties for the masses are listed in Tables 5 and 6. For the cross-section and cross- 
section times branching fraction measurements, the smallest allowed relative precision is 
raised to 2 % for the fit including theoretical uncertainties (as opposed to 1 % for the fit 
without theoretical uncertainties). The full list of observables used in the fit and their 
uncertainties can be obtained from [24]. 



4 Fit Results 

The input observables described in Section 3.2 are used to determine the SUSY Lagrangian 
parameters in a global fit under the assumptions mentioned in Section 3.1. In total 18 SUSY 
parameters remain to be fitted. In addition to those, the top mass mt is fitted, since it has 
a relatively large uncertainty and strongly influences parts of the MSSM observables. Thus 
19 Lagrangian parameters are simultaneously determined in this fit. 

As shown in Table 1 all parameters are perfectly reconstructed at their input values. Due 
to the fact that the input observables are unsmeared in this fit, the final x 2 is close to zero 
at x 2 = 2-1 x 1(T 5 . 
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Parameter 


"True" value 


Fit value 


Uncertainty 


Uncertainty 








(exp.) 


(exp.+theor.) 


tan/5 


10.00 


10.00 


0.11 


0.15 




400.4 GeV 


400.4 GeV 


1.2 GeV 


1.3 GeV 




-4449. GeV 


-4449. GeV 


20. GeV 


29. GeV 


Me, 


115.60 GeV 


115.60 GeV 


0.13 GeV 


0.43 GeV 


Mf R 


109.89 GeV 


109.89 GeV 


0.32 GeV 


0.56 GeV 




181.30 GeV 


181.30 GeV 


0.06 GeV 


0.09 GeV 


M fL 


179.54 GeV 


179.54 GeV 


0.12 GeV 


0.17 GeV 


x t 


-565.7 GeV 


-565.7 GeV 


6.3 GeV 


15.8 GeV 


x h 


-4935. GeV 


-4935. GeV 


1207. GeV 


1713. GeV 




503. GeV 


504. GeV 


12. GeV 


16. GeV 


M ln 


497. GeV 


497. GeV 


8. GeV 


16. GeV 


M iR 


380.9 GeV 


380.9 GeV 


2.5 GeV 


3.7 GeV 


M« L 


523. GeV 


523. GeV 


3.2 GeV 


4.3 GeV 


M iL 


467.7 GeV 


467.7 GeV 


3.1 GeV 


5.1 GeV 


Mi 


103.27 GeV 


103.27 GeV 


0.06 GeV 


0.14 GeV 


M 2 


193.45 GeV 


193.45 GeV 


0.08 GeV 


0.13 GeV 


M 3 


569. GeV 


569. GeV 


7. GeV 


7.4 GeV 


™A run 


312.0 GeV 


311.9 GeV 


4.3 GeV 


6.5 GeV 


m t 


178.00 GeV 


178.00 GeV 


0.05 GeV 


0.12 GeV 


Correspondinj 


y values for the trilinear coupling 


;s: 




A T 


-445. GeV 


-445. GeV 


40. GeV 


52. GeV 


A t 


-526. GeV 


-526. GeV 


6. GeV 


16. GeV 


A h 


-931. GeV 


-931. GeV 


1184. GeV 


1676. GeV 



X 2 for unsmeared observables: 2.1 x 10 



Table 1: The Fittino fit result for the SPSla' inspired scenario. The left column shows 
the values predicted by SPheno version 2.2.2 for this scenario, the second column 
exhibits the central values of the fit with unsmeared input observables, the third 
column displays the parameter uncertainties for the fit with experimental uncer- 
tainties only. The last column shows the parameter uncertainties for the fit with 
experimental and theoretical uncertainties. In addition to the fitted parameters, 
the corresponding values for the trilinear couplings are given as obtained from 
reversing the parameter transformations (6)-(8), taking correlations into account. 

4.1 Extraction of the Fit Uncertainties and Correlations 

After a successful convergence of the simulated annealing algorithm to the input parameter 
values, the fit uncertainties are evaluated by carrying out many individual fits with input 
values randomly smeared within their uncertainty range using a Gaussian probability density. 
The complete covariance matrix and the correlation matrix are derived from the spread of 
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Parameter 


(RMS - a)/ a 


Parameter 


(RMS - a)/(7 


tan/5 


0.01 


M iR 


0.07 




0.06 


0.02 




0.03 


M- qL 


0.01 




0.02 


M iL 


0.04 


M fR 


0.00 


Mi 


0.01 




0.02 


M 2 


-0.01 


M fL 


0.01 


M 3 


0.04 


x t 


0.03 


m\ 

^run 


0.14 


x h 


0.17 


m t 


0.01 




0.04 







Table 2: Comparison of the RMS and the Gaussian width a fitted to the parameter distri- 
butions. 

the fitted parameter values. 

For the case without and for the case with theoretical uncertainties about 1000 fits are 
performed. For each of those fits, the \ 2 is minimized. For a large and complex parameter 
space with large correlations among the parameters this method has turned out to be more 
robust than a MINOS error analysis. 

An example of the outcome of this procedure for the fit (with experimental uncertainties 
only) from Section 4 is shown in Figure 2 for the parameters tan/5, Mi, ha po1c and Xb for 
the ~1000 independent fits. All distributions apart from ?tia po1c and X h agree well with 
Gaussians. 

The difference of the RMS of the parameter distribution and the width of a Gaussian 
fitted to the distribution (normalized to the Gaussian width) is shown in Tab. 2 for the 19 
parameters. The distributions showing the largest disagreement are ?tia po1c and At,, indicating 
that a parabolic error assumption for these parameters is not completely correct. 

The parameter uncertainties extracted from the fit value distributions are shown in 
Tab. 1, while the corresponding correlation matrix is displayed in Tab. 7 and 8. Addi- 
tionally, a graphical representation of the relative uncertainties of the individual parameters, 
both with and without theoretical uncertainties, is given in Figure 1. 

Figure 3 shows the distribution of the x 2 values for these ~1000 independent fits for the 
case without theoretical uncertainties. The mean \ 2 obtained from a fit of the x 2 distribution 
to the observed distribution is 128.0, in agreement with the expectation of 129.0 ± 0.7. This 
shows that the fits converge well to the true minimum of the \ 2 surface for the smeared 
observables, implying that the uncertainties extracted from the toy fit value distributions 
are correct. 
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Figure 2: Examples of the toy fit value distributions for ~1000 independent fits with ob- 
servables smeared within their uncertainties. 
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Figure 3: % 2 distribution for ~1000 independent fits with observables smeared within their 
uncertainties. The mean \ 2 from a x 2 function fitted to the observed distribution 
of 128.0 agrees with the expectation of 129.0 ± 0.7. 
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more precise calculations would be beneficial. 

The least constrained variables are the mixing parameters. Especially has only small 
influence on the Higgs sector and the sbottom masses, resulting in a relative precision of 
only 25 % (corresponding to 130 % relative uncertainty for A\>). The situation is better for 
X T and X t . The high precision in the chargino and slepton sector as well as the occurrence 
of the decay channel H + — * v Th f^ allows to determine X T to an accurary of less than 0.5 % 
(corresponds to 9 % for A T ). The achieved relative precision of 1 % for X t (1 % for A t ) 
is possible due to the large mass splitting in the stop sector and its influence on the Higgs 
sector. Searching for sensitive observables to constrain the mixing parameters to a larger 
extent is important to improve their precisions. 

The benefit of a global fit is that the full correlation matrix of all parameters is available 
(see Appendix B). As expected correlations within sectors can be strong. The parameters of 
the gaugino sector tan ft, [i, Mi and M 2 , for example, reveal correlations up to 69 %. But even 
between different sectors the correlations cannot be neglected. One example is M^ R from the 
squark sector whose correlation coefficient with M gfl from slepton sector amounts to —0.79. 
Even the relatively robust gaugino sector is affected by such inter-sector correlations. M 2 , 
for example, has a correlation of 20 % with M^ L from the sfermion sector. This underlines 
the importance of fitting all SUSY parameters simultaneously. Neglecting these inter-sector 
correlations causes incorrect fit uncertainties in fits restricted to specific sectors. 

4.3 Fits in Subspaces 

This section is devoted to a check of the influence of restricted fits on the precision of 
the reconstructed parameters and their uncertainties. It might be interesting to study the 
possibility to extract limited information on the SUSY parameters in a certain subspace 
only, because the input described in Section 3.2 will not be completely available at the 
beginning of ILC running at lower energies. Therefore fits with inputs restricted to the low- 
energy running phase of ILC have been performed. But this entails fixing parameters which 
cannot be determined with the low-energy data. The following excerpt of observables from 
Section 3.2 has been used for the fit: 

• All SM observables. 

• All mass measurements in the gaugino and slepton sector from ILC running at y/s = 
500 GeV. 

• All chargino, neutralino and e cross-section times branching fraction measurements at 
y/s = 400 GeV and ^1 = 500 GeV. 

No Higgs sector parameters and observables are used in the fit, since they have large depen- 
dencies on squark sector parameters such as X t and Mj. 
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Parameter 


"True" value 


Fit with correctly 


Fit with incorrectly 


Uncertainty 






fixed parameters 


fixed parameters 




Fixed parameters 


x t 


-565.7 GeV 


-565.7 GeV 


-30.0 GeV 


fixed 


x h 


-4934.8 GeV 


-4934.8 GeV 


-4000.0 GeV 


fixed 




501.6 GeV 


501.6 GeV 


600.0 GeV 


fixed 


M ~ bR 


497.4 GeV 


497.4 GeV 


600.0 GeV 


fixed 




503.9 GeV 


503.9 GeV 


500.0 GeV 


fixed 




523.2 GeV 


523.2 GeV 


500.0 GeV 


fixed 




467.7 GeV 


467.7 GeV 


500.0 GeV 


fixed 


M 3 


568.9 GeV 


568.9 GeV 


700.0 GeV 


fixed 


™A run 


312.0 GeV 


312.0 GeV 


400.0 GeV 


fixed 




178.0 GeV 


178.0 GeV 


178.0 GeV 


fixed 


Fitted parameters 


tan p 


10.00 


10.00 


11.1 


0.47 




400.39 GeV 


400.388 GeV 


388.3 GeV 


3.1 GeV 


x T 


-4449.2 GeV 


-4449.2 GeV 


-4447.8 GeV 


37.2 GeV 


Me R 


115.60 GeV 


115.602 GeV 


113.74 GeV 


0.06 GeV 


M fH 


109.89 GeV 


109.89 GeV 


107.77 GeV 


0.48 GeV 


M, L 


181.30 GeV 


181.304 GeV 


181.76 GeV 


0.04 GeV 


M fL 


179.54 GeV 


179.54 GeV 


179.99 GeV 


0.14 GeV 


M 1 


103.271 GeV 


103.271 GeV 


103.11 GeV 


0.05 GeV 


M 2 


193.446 GeV 


193.445 GeV 


193.49 GeV 


0.12 GeV 


x 2 




1.8 x 10~ 5 


5.89 





Table 3: Fittino fits of a parameter subspace using inputs from the low-energy running of 
ILC only. 



The fixed and fitted parameters are listed in Table 3. All fixed parameters are set to 
the correct values in a first fit. As shown in Table 3 the central values of the fit reproduce 
the correct parameter values. In a second step the fit is performed with all fixed parameters 
set to modified values. The chosen numbers are estimates in a situation where the strongly 
interacting sector has not been precisely measured yet. The quality of the fit results for this 
second fit depends on the chosen set of fixed and fitted parameters. 

This test shows that subspace fits can provide results of the correct order of magnitude. 
Subspace fits are justified at the beginning of the exploration of SUSY phenomena where 
experimental uncertainties are large and only limited data is available for SUSY parameter 
analysis. Depending on the chosen set of fixed and fitted parameters, subspace fits might 
yield parameter values whose biases are several times as large as their uncertainties obtained 
from the fit. This unfortunate feature is a consequence of neglected correlations. The errors 
obtained from subsector fits do not express a probability to find the true parameter values 
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within the uncertainty range. Uncovering these shifts by looking at the x 2 of the fit is 
unlikely in reality since a degradation of A% 2 = 5.89, which is the largest x 2 enhancement 
observed in above fits, might still yield a reasonable x 2 value for the fit with its 85 degrees 
of freedom. In any case, as this example shows, a precision determination of the SUSY 
Lagrangian parameters requires a comprehensive global fit using inputs from all sectors of 
the theory and fitting parameters from all sectors of the theory. 



5 Extrapolation to High Energy Scales 



As soon as the complete set of parameters is known at the electroweak scale, they can 
be extrapolated to high energy scales to test ideas concerning SUSY breaking or grand 
unification within the bottom- up approach [7]. In this way one exploits the experimental 
information to the maximum extent possible without any assumptions on the structure of 
the high scale theory except that there are no new particles between the electroweak scale 
and the high scale 1 . 



a) Mf 1 [GeV- 1 
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Figure 4: Running of a) the gaugino mass parameters, b) first generation scalar mass 
parameters and Mfj 2 and c) third generation scalar mass parameters and M\ 
in SPSla'. Full bands: only experimental errors are taken into account; dashed 
lines: today's theoretical errors are taken into account as a conservative estimate. 



The parameters at the high energy are related to the electroweak scale parameters by 
renormalization group equations (RGEs). Here we use two-loop RGEs for the evolution [23]. 
The evolution of the parameters from the electroweak scale to the GUT scale are shown 
in Fig. 4 where the GUT scale is defined as the scale where the SU{2) gauge coupling g 2 
coincides with U(l) gauge coupling gi (using the proper normalisation gi = ^5/3g') and 

■"■In the case that additional heavy particles, e. g. right handed neutrinos, one is even able to get information 
on their mass scale in this approach as has been shown in [34]. 
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Figure 5: Running of third generation scalar mass parameters and M] Ii in SPSla' assuming 
that A\> is known with an 'accuracy' of 50%. 

the value for the electroweak scale is given by y/m^m^. The latter value is motivated by the 
expectation that at this scale the effect of missing higher order corrections in the calculation 
of mh, which is most probably the most precisely measured mass, is minimized. 

In Fig. 4 a the evolution for the gaugino parameters Mf is presented. They are clearly 
under excellent control allowing for precise tests concerning their unification at Mqut- Note, 
that their 'unification point' is in principle independent of the one for the gauge couplings 
allowing for a cross check if both sets of parameters indeed meet at the same point. 

In Fig. 4 b and c the running of the scalar mass parameters are shown. The slepton mass 
parameters are under excellent control whereas the squark and the Higgs mass parameters 
are signifcantly worse. Part of this difference can be traced back to the structure of the 
RGEs as explained in detail in the second paper of [7]. In addition the analysis used here 
does not coincide with ones presented in [7] and [35]: In [7] and [35] it has been assumed 
that all A parameters meet within 1-sigma as up to now no good proposal exists for a precise 
measurement of A^. As we did not use this assumption in this study we find significantly 
larger errors for the third generation squark mass parameters and the Higgs parameters. 
Using this assumption or, equivalently, having tools at hand which allow for a better deter- 
mination of Ab, the error bands get significantly reduced as can be seen in Fig. 5 where we 
show the running of the same parameters as in Fig. 4 c assuming that A^ is known within 
50 % at the electroweak scale. 

In Fig. 4 we also show the case that today's theoretical uncertainties are taken into 
account (dashed lines). In case of the gaugino mass parameters the effect is so small that 
it cannot be seen in the plot (see also Tab. 1). Note that the smallness of the theoretical 
uncertainty on the gluino mass is due to the fact that the squarks have a similar mass. It gets 
larger if there is a significant hierarchy between these particles, e. g. in focus point scenarios. 
One sees that the main effect is on the parameters of the third generation sfermions and 
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SPSla' value 


Fitted value 


Alhc+ilc 


AlHC only 


tan (3 


10.000 


10.000 


0.036 


1.3 


M (GeV) 


70.000 


70.000 


0.070 


1.4 


M 1/2 (GeV) 


250.000 


250.000 


0.065 


1.0 


A) (GeV) 


-300.0 


-300.0 


2.5 


16.6 



Table 4: Fit results for fits within mSUGRA scenario. The meanings of the columns are 
(starting from the left): SPSla' values, the fitted mSUGRA parameters, parameter 
uncertainties for a fit to LHC+ILC observables and parameter uncertainties for 
a fit to "LHC only" observables. In both cases theoretical uncertainties are not 
included. 

the Higgs mass parameters. The main reason for this is that the error on and M^ r get 
considerably larger once the theoretical uncertainties are taken into account. This figure 
clearly shows that there is a strong need to further increase accuracies of the theoretical 
predictions. 

6 Fit Within mSUGRA Framework 

Instead of determining the full spectrum of the general low-energy Lagrangian parameters 
and inferring the high-scale mechanisms from them, the parameters of a more restrictive high- 
scale scenario can be fitted directly to the observables. This section presents an example for 
the mSUGRA scenario SPSla', using the same set of observables as described in Section 3.2. 
The parameters of the scenario are tan/3, M , M1/2, A an d sign(/x). The latter has been 
fixed to 1 for this exercise. 

Since there are only four continuous parameters, the fit converges without complex fitting 
techniques and the uncertainty determination is stable. The fit results for the simulated 
LHC+ILC measurements from Section 4 and observables set to the exact prediction of 
the mSUGRA model are given in Table 4, along with results for a fit just using the LHC 
observables in the fit. As an additional LHC measurement, the ratio of Higgs couplings as 
simulated in [36] has been used. For the LHC+ILC measurements, due to the much more 
constrained theoretical parameter space, the precision on the parameters is much better than 
for the general low-energy model (see Table 1). For the same reason, the LHC alone is able 
to perform well constrained parameter measurements with uncertainties typically one to two 
orders of magnitude larger than in the case of the combination of both. 

An important test is, however, to determine how well the experimental results can be 
used to determine deviations from a specific high-scale model such as mSUGRA. Therefore 
the mSUGRA high-energy parameters have been fitted to the observables from the SPSla' 
inspired scenario used in Section 3, which is only slightly modified to obtain exact unification 
in the first two generations of squarks and sleptions. For this case, the total x 2 s from the 
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two fits indicate that the LHC+ILC measurements are precise enough to discriminate the 
modified scenario from exact mSUGRA (% 2 = 3600, 144 n.d.f.), while the LHC measurements 
alone do not allow this insight (x 2 = 0.5, 25 n.d.f.). 

7 Conclusions 

The Lagrangian parameters of the MSSM assuming universality for the first and second 
generation and real parameters but without assumptions on the SUSY breaking mechanism 
are correctly reconstructed without usage of a priori information. This has been achieved 
using precision measurements at the LHC and ILC as input to a global fit exploiting the 
techniques implemented in the program Fittino. 

Most of the Lagrangian parameters can be determined to a precision around the percent 
level. For some parameters an accuracy of better than 1 per mil is achievable if theoretical 
uncertainties are neglected. As a (very) conservative estimate we have re-done the fit using 
today's theoretical uncertainties. As expected it turns out that they can significantly dete- 
riorate the precision of the Lagrangian parameter determination. More work is needed to 
improve the accurary of theoretical predictions in order to fully benefit from the experimen- 
tal precision. The only parameter which cannot be strongly constrained by the observables 
used in the presented fit is resulting in a precision of only 25 % for this parameter which 
translates into an uncertainty of 130 % on A^. 

The strong (inter-sector) correlations between the fitted parameters underline the im- 
portance of performing a global fit including all sectors of the theory. As shown, fits in 
parameter subspaces can lead to significant biases and wrong uncertainty estimates. 

The fitted parameters at the electroweak scale have been extrapolated to high scales 
using the bottom-up approach taking into account all correlations between the errors. We 
have found that the gaugino sector as well as the slepton sector can be reconstructed rather 
precisely allowing for stringent tests of the underlying model in these sectors. The accuracy 
deteriorates considerably for squark mass parameters and the Higgs mass parameters. The 
results clearly show how far the squark mass parameters can deviate from the universality 
assumption indicated by the measurements in the slepton sector. We have also performed a 
top-down fit of the high-energy mSUGRA parameters to combined LHC+ILC observables. 
Here we have found that they are very sensitive to even tiny deviations from the mSUGRA 
hypothesis. The discriminative power of "LHC only" inputs is significantly weaker. 
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A List of Input Observables to the SPSla' Inspired Fit 



Observable Value Exp. uncertainty 

"m^ 80.3371 GeV 0.039 GeV 

m z 91.1187 GeV 0.0021 GeV 

m t 178.0 GeV 0.05 GeV 

m h 112.888 GeV 0.05 GeV 

m Apolc 374.228 GeV 1.3 GeV 

m H P ° 374.464 GeV 1.3 GeV 

m H ± 383.131 GeV 1.1 GeV 

m- qL 561.539 GeV 9.8 GeV 

m- qR 543.35 GeV 11.0 GeV 

m~ bl 502.059 GeV 5.7 GeV 

m- b2 541.81 GeV 6.2 GeV 

ml 365.819 GeV 2.0 GeV 

m e - L 190.209 GeV 0.2 GeV 

124.883 GeV 0.05 GeV 

190.237 GeV 0.5 GeV 

m~ m 124.837 GeV 0.2 GeV 

m fl 107.292 GeV 0.3 GeV 

m f2 195.290 GeV 1.1 GeV 

m g 603.639 GeV 6.4 GeV 

m^o 97.7662 GeV 0.05 GeV 

X.2 



m e R 



184.345 GeV 0.08 GeV 

7tv> 404.134 GeV 4.0 GeV 

A3 
A4 



417.037 GeV 2.3 GeV 

m-± 184.132 GeV 0.55 GeV 

m~± 418.495 GeV 3.0 GeV 

Edge 3 with m^o, m & , m^o 449.679 GeV 4.9 GeV 

Edge 3 with m^ R , rriq L , m^o 390.285 GeV 3.35 GeV 

Edge 4 with m^o, m^o, m^ R , rriq L 329.831 GeV 4.2 GeV 

Edge 5 with m^o, m x9,i m jj.R-> m <iL 218.529 GeV 3.44 GeV 

Table 5: Simulated measurements at LHC and at ILC running at 400 GeV, 500 GeV and 
1 TeV center-of-mass energy. The values of the observables are taken from the 
prediction of SPheno version 2.2.2 for the SPSla' inspired scenario. The experi- 
mental uncertainties on the masses are taken from [15]. The definition of the edge 
numbers can be found in [9]. 
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Observable Theor. uncertainty 

m t 

m H 

m H ± 

m Q L 
m Q R 
m ' bl 
m l 2 



m iR 



m- g 

X-l 
A3 

m xf 

Edge 3 with m^o, rriq L , m^o 
Edge 3 with rrip jR , rrig L , m^o 
Edge 4 with m^o, m^o, m^ R , rriq L 
Edge 5 with m^o, m^o, m g - L 
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.5 


GeV 





.5 


GeV 


1 


4 


GeV 





.4 


GeV 


1 


.2 


GeV 


1 


.2 


GeV 


1 


.2 


GeV 


1 


.3 


GeV 


1 


.3 


GeV 


4 


.5 


GeV 


3 


.9 


GeV 


3 


.3 


GeV 


2 


49 GeV 



Table 6: Theoretical uncertainties on the masses (taken from [17]) and the induced ones on 
the edge variables. 
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B Correlation Matrix of the SPSla' Inspired Fit 

Table 7 and 8 show the full correlation matrix from the fit for the SPSla' inspired scenario 
described in Section 3. 



Parameter 


tan /3 




X T 






fn 


M s, 


ft f 

Mf T 


x t 


x b 




tan P 


1 


.000 


—0.509 


—0.238 


0.042 





032 


-0 


125 


—0.321 





030 


-0 


569 


0.033 


A 4 


— u 


ouy 


1.000 


0.014 


—0.038 


— 


028 


n 
U 


oUo 


0.308 


— U 




n 
U 


101 


0.031 


Y 


-0 


,238 


0.014 


1.000 


0.055 


_Q 


431 


-0 


,076 


0.225 


-0 


137 


-0 


091 


—0.058 


m Sr 





.042 


-0.038 


0.055 


1.000 





.369 


-0 


631 


-0.329 





026 


-0 


,109 


-0.793 


M fR 





.032 


-0.028 


-0.431 


0.369 


1 


.000 


-0 


,254 


-0.339 





062 





,074 


-0.362 




-0 


.125 


0.303 


-0.076 


-0.631 


-0 


.254 


1 


,000 


0.368 


-0 


,048 





,112 


0.557 


M fL 


-0 


321 


0.308 


0.225 


-0.329 


-0 


.339 





,368 


1.000 


-0 


,023 





121 


0.320 


x t 





.030 


-0.064 


-0.137 


0.026 





.062 


-0 


,048 


-0.023 


1 


000 


-0 


360 


0.097 


x h 


-0 


.569 


0.151 


-0.091 


-0.109 





.074 





,112 


0.121 


-0 


,360 


1, 


,000 


-0.090 


m 9r 





033 


0.031 


-0.058 


-0.793 


-0 


.362 





,557 


0.320 





,097 


-0 


,090 


1.000 


M l R 


-0 


.049 


0.037 


-0.016 


0.193 





.080 


-0 


133 


-0.011 





,076 





004 


0.091 


M t R 

m 9l 


-0 


166 


0.125 


-0.011 


-0.335 


-0 


141 





,327 


0.210 


-0 


371 





307 


0.082 


-0 


.065 


0.025 


-0.035 


0.192 





.069 


-0 


,084 


-0.001 





144 


-0 


,014 


0.076 


Mt L 


-0 


.005 


-0.094 


0.054 


0.109 





.001 


-0 


,138 


-0.066 


-0 


,270 





092 


0.070 


Mi 


-0 


.147 


-0.263 


0.118 


0.169 





.051 


-0 


,177 


0.032 





,044 





130 


-0.244 


M 2 





.229 


-0.694 


-0.299 


-0.047 





.062 





,087 


0.024 





,128 





,042 


-0.058 


Ah 





.044 


0.003 


-0.036 


0.078 





.065 


-0 


,046 


-0.061 





,076 


-0 


,070 


-0.180 


m A run 


-0 


.594 


0.143 


-0.048 


-0.108 





.056 





,110 


0.130 


-0 


,223 





,972 


-0.076 


m t 


-0 


.004 


-0.013 


0.004 


-0.005 





.018 


-0 


,018 


0.039 


0, 


,186 


-0 


,029 


0.011 



Table 7: Correlation matrix of the Fittino SPSla' motivated fit, part I. 



Parameter 


Mr 
bn 


M f 


m 9l 


M f 


Mi 


M 2 


M 3 


m A run 


mt 


tan /3 


-0.049 


-0 


166 


-0 


,065 


-0.005 


-0 


,147 





229 


0.044 


-0 


,594 


-0.004 


A 4 


0.037 





125 





,025 


-0.094 


-0 


,263 


-0 


,694 


0.003 





,143 


-0.013 


x T 


-0.016 


-0 


Oil 


-0 


,035 


0.054 





,118 


-0 


,299 


-0.036 


-0 


,048 


0.004 


m Sr 


0.193 


-0 


,335 





,192 


0.109 





,169 


-0 


,047 


0.078 


-0 


,108 


-0.005 


Mf R 


0.080 


-0 


141 





,069 


0.001 





,051 





062 


0.065 





,056 


0.018 


Me L 


-0.133 





,327 


-0 


,084 


-0.138 


-0 


,177 





,087 


-0.046 





,110 


-0.018 


Mf L 


-0.011 





,210 


-0 


,001 


-0.066 





,032 





,024 


-0.061 





130 


0.039 


x t 


0.076 


-0 


,371 





,144 


-0.270 





044 





,128 


0.076 


-0 


,223 


0.186 


x h 


0.004 





,307 


-0 


,014 


0.092 





130 





,042 


-0.070 





,972 


-0.029 


Mq R 


0.091 





,082 





,076 


0.070 


-0 


,244 


-0 


,058 


-0.180 


-0 


,076 


0.011 


M l R 


1.000 





,009 


-0 


,017 


-0.077 


-0 


,034 





,024 


-0.105 





,026 


0.012 


M l R 


0.009 


1 


,000 


-0 


036 


-0.428 


-0 


,048 





,045 


-0.236 





,249 


-0.063 


M~ qL 


-0.017 


-0 


,036 


1 


,000 


0.057 





,006 


-0 


,200 


-0.275 





Oil 


0.017 


Mr 

Mi 


-0.077 


-0 


,428 





,057 


1.000 


-0 


002 


-0 


,140 


-0.225 





,106 


-0.081 


-0.034 


-0 


,048 





,006 


-0.002 


1 


,000 





293 


0.055 





,141 


-0.060 


M; 


0.024 





,045 


-0 


,200 


-0.140 





,293 


1, 


,000 


0.109 





,040 


-0.005 


M 3 


-0.105 


-0 


,236 


-0 


,275 


-0.225 





,055 





,109 


1.000 


-0 


,086 


0.008 


m A run 


0.026 





,249 





011 


0.106 





141 





040 


-0.086 


1 


,000 


-0.008 


mt 


0.012 


-0 


,063 





,017 


-0.081 


-0, 


,060 


-0 


,005 


0.008 


-0 


,008 


1.000 



Table 8: Correlation matrix of the Fittino SPSla' motivated fit, part II. 
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